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Abstract — We derive in a direct and rather straightforward 
way the null controllability of a 2-D heat equation with 
boundary control. We use the so-called flatness approach, which 
consists in parameterizing the solution and the control by the 
derivatives of a "flat output". This provides an explicit control 
law achieving the exact steering to zero. Numerical experiments 
demonstrate the relevance of the approach. 



I. Introduction 

The controllability of the heat equation was first con- 
sidered in the 1-D case, [l]-[3]), and very precise results 
were obtained by the classical moment approach. Next 
using Carleman estimates and duality arguments the null 
controllability was proved in [4], [5] for any bounded domain 
in R^, any control time T, and any control region. This 
Carleman approach proves very efficient also with semilinear 
parabolic equations, [4]. 

By contrast the numerical control of the heat equation (or 
of parabolic equations) is in its early stage, see e.g. [6]- 
[8]. A natural candidate for the control input is the control 
of minimal L^— norm, which may be obtained as a trace 
of the solution of the (backward) adjoint problem whose 
terminal state is the minimizer of a suitable quadratic cost. 
Unfortunately its computation is a hard task [8]; indeed the 
terminal state of the adjoint problem associated with some 
regular initial state of the control problem may be highly 
irregular, which leads to severe troubles in the numerical 
computation of the control function. 

All the above results rely on some observability inequal- 
ities for the adjoint system. A direct approach which does 
not involve the adjoint problem was proposed in [2], [3], 
[9]. In [2] a fundamental solution for the heat equation 
with compact support in time was introduced and used to 
prove null controllability. The results in [2], [10] can be 
used to derive control results on a bounded interval with 
one boundary control in some Gevrey class. An extension 
of those results to the semilinear heat equation in ID was 
obtained in [9] in a more explicit way through the resolution 
of an ill-posed problem with data of Gevrey order 2 in t. 

In this paper we derive in a straightforward way the null 
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controllability of the 2-D heat equation 

Ot{t, Xi,X2) = A0{t,Xi,X2) (1) 

0:,,{t,0,X2)=0 (2) 

0,,{t,L,X2) = (3) 

0,,{t,xi,0) = (4) 

Ox^{t, xi,l) = u{t,xi) (5) 

with initial condition 

0{0,xi,X2) = 6o{xi,X2); 

here (xi,X2) belongs to the rectangle Q := (0,L) x (0, 1), 
t e (0,T), and Oq G L'^{Q). This system describes the 
dynamics of the temperature 6 in an insulated metal strip 
where the control u is the heat flux at one end. 

More precisely given any final time T > and any initial 
state Oq e L'^{ft) we provide an explicit very regular control 
input u such that the state reached at time T is zero, i.e. 

0{T, xi,X2) =0, xen. 

We use the so-called ^am^^^ approach, [11], which consists 
in parameterizing the solution and the control u by the 
derivatives of a "flat output" y (section |Il]); this notion 
was initially introduced for finite-dimensional (nonlinear) 
systems, and later extended to (in particular) parabolic PDEs, 
[12]-[15]. Choosing a suitable trajectory for this flat output y 
then yields an explicit series for a control achieving the 



exact steering to zero (section III). Numerical experiments 



demonstrate the relevance of the approach (section [rv). 

This paper extends [16] to the 2-D case, with moreover a 
more elegant and efficient construction of the control. This is 
probably the first example using flatness for motion planning 
of a "truly" 2-D PDE. 

In the sequel we will consider series with infinitely many 
derivatives of some functions. The notion of Gevrey order 
is a way of estimating the growth of these derivatives: we 
say that a function y G C^([0, T]) is Gevrey of order s >0 
on [0, T] if there exist positive constants M, R such that 



y^'Ht) 



< M 
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RP 



yte [o,T], Vp>o. 



More generally if K C M^ is a compact set and y is 3. 
function of class C^ on K (i.e. y is the restriction to K 
of a function of class C^ on some open neighbourhood ft 
of K), we say y is Gevrey of order si in xi, S2 in X2,. . .,sn 
in xjsf on K if there exist positive constants M, i?i, ..., Rn 
such that for all x e K and p G N^ 
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By definition, a Gevrey function of order s is also of order 
r for r > 5. Gevrey functions of order 1 are analytic (entire 
if 5 < 1). Gevrey functions of order s > 1 have a divergent 
Taylor expansion; the larger s, the "more divergent" the 
Taylor expansion. Important properties of analytic functions 
generalize to Gevrey functions of order s > 1: the scaling, 
addition, multiplication and derivation of Gevrey functions 
of order s > 1 is of order s, see [17], [18]. But contrary 
to analytic functions, functions of order s > 1 may be 
constant on an open set without being constant everywhere. 
For example the "step function" 



where the a^'s are functions yet to define. Plugging this 
formal solution into ([T]) we find 



2],^[^^+2(^^^i) - i^t -dljai{t,xi)] =0, 



Mt) '-- 



-ii-ty 



K e 



-(i-t)- 



if t < 
ift>l 

if t g]0,1[ 



where k = {s — 1) ^, is Gevrey of order s on [0, 1] (and 
in fact on R); notice ^s{0) = 1, (/>s(l) = and (/)i'^(0) = 
#(1) =0 for alH> 1. 

In conjunction with growth estimates we will repeatedly 
use Stirling's formula n! ^ {n/e)^V27rn. 

II. The heat equation is "flat" 

Let {ej)j>o be an orthonormal basis of L^(0,L) such 
that each function Cj is an eigenfunction for the Neumann 
Laplacian on (0,L), ie is a solution of 

-e-{xi) = Xjej{xi), xi e (0,L) 
e;.(0)=0 
e;.(L)=0. 
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For instance A. 



and 



eo(xi) = L~^ 
Decompose then 6>(t,xi,0) as 



2 fJ7rxi\ 



, J>1. 



i9(t, xi,0) = ^Zj{t)ej{xi). 

j>0 

We claim ^-^ is "flat" with z{t) := {zj{t)) .^^ as a 
flat output, which means there is (in appropriate spaces of 
smooth functions) a 1 — 1 correspondence between arbitrary 
functions t \-^ z{t) and solutions of ([T])-([5]). 

On the one hand a solution of ([l])-([5]) uniquely defines 
the Zfc's hence z by 

pL pL 

I 0{t^xi^O)ek{xi)dxi = 2_.^j{t) / ej{xi)ek{xi)dxi 
^.>o Jo 



I 



On the other hand let us seek a formal solution in the form 

0{t, xi, X2) = ^ -^ai{t, xi), 



i>0 



i>0 

hence a^+2 = {dt — d^Jai for all i > 0. Moreover 
ao(t,xi) = 6>(t,xi,0) and by ^ ai{t^xi) = 0. Therefore 
for aU i > 

a2i{t,xi) = {dt -dl^yao 

= Y.(9t-dlyz,{t)e,{xr) 

j>0 

= ^e,{x^){dt^X,yz,{t) 

j>0 

= Y.^,{x,)e-^^'yf{t), 

j>0 

where we have set yj{t) := e^^^Zj{t). Clearly 



^2i 



^(t,xi,X2) = ^e-^^S-(^i)E^fW(^ 



j>0 



i>0 



4t,xo = E^"'^S(^i)E7^-^^fw- 



(6) 

(7) 



is then a formal solution of ([T])-([5]) uniquely defined by 
the ?/j's hence by z. 

We now give a precise meaning to this formal solution by 
restricting the ?/j's to be Gevrey of order s G [1, 2). 

Theorem 1: Let s G [1,2), < ti < ^2 < 00, and a 
sequence y = {yj)j>o in C^([0,T]) satisfying for some 
constants M,R> 



R^ 



yfi^) <^^' vz,j>o, yte[h,t2]. 



Then the formal solution defined by ([6]) is Gevrey of order s 

in t, 1/2 in xi, and s/2 in X2 on [ti, ^2] xl]; as a consequence 
the formal control u defined by ^ is also Gevrey of order s 
in t and 1/2 in xi on [^1,^2] x [0,L]. 

Proof: We must prove the formal series 

2i>n j>0 



where 



2i-n 

^'^■•" (2i-n)! 



{-X,ye,{x,)dr{e-'^'yf{t)). 



are uniformly convergent on [^1,^2] x ft with growth esti- 
mates of the form 



\drdi[d:^j{t,x,,x2)\<c 



m\^ (2/)! 2 n!2 
^ Rf ~B^ 



for some positive constants i?o , ^1 , ^2 • In what follows C 
denotes a generic constant independent of m, /, n, t,Xi,X2 



that may vary from line to line. By assumption 
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^.(-i)Eu br^w^r-'^-'^* 



(2i-n)! •'' " ^\k 

^ ' k=0 

ml 
(2i - n)! ^ k\(i^ k)\ Ri+^ "^ 






m! (i + fc)!* 



A'- 



l-\-m — k—\jti 



Let to G (0,ti) and^ := ti-to. Using (p+g)! < 2^+^^!^! 
and x^e~^ < p! for x > 



k\R^+^ - k\R^+^ 



"-'i^y 



,,.„-V-M<<4S=i^ <,!(„. -«!(H)'^'"-', 



t' 



so we have 



l^i.-l <e-^^''M 



'''{^) ,. ../2\^+-^/2M\^ 



/^x L+m ^^ / 



L (2i - n) 
Now by Stirling's formula 



«0~<-)'*^^-e|i, 



fe=0 



2R 



where i?i < V^. Likewise using also (p + g)! < 2^+^^!^^! 



'•'{i) fry (2i)!^ (^i)^ 



{2i-n)\ \RJ {2i-n)i\{2i-n)\^-i2^ 
'2^\'^ /n7T\i / /2i — n\i 



^y (t) ^+(^ 



n!2 



{2i-n)\^-2 



<C 



n'2 -^ ^ \^ n y 

^ (2i-n)!i-t V i^ 



A2i-. 



IIL Null controllability 

In this section we derive an explicit control steering the 
system from any initial state 6o G L'^{ft) at time to the 
final state at time T > 0. Two ideas are involved: on 
the one hand thanks to the flatness property it is easy, by 
selecting a suitable trajectory for the flat output, to find a 
control achieving the steering to zero starting from a certain 
set of initial conditions; on the other hand thanks to the 
regularizing property of the heat equation this set is reached 
from any Oq G L'^{Q) when applying first a zero control for 
some time (lemma [T]). 

Lemma 1: Let 6q G L'^{Q) and r > 0. Consider the final 
state 0r{xi^X2) := ^(t, Xi,X2) reached when applying the 
control u{t) := 0, t G [0, r], starting from the initial state Oq. 

Then 0^ can be expanded as 



0r{xi,X2) = ^e ^'^ej{xi)^yj. 



y,2i 



i>o 



i>0 



(2i)!' 



with 



\y,A<c{i+') 



rJ T 



where R2 < y ^\ notice ^2i>n ^^i-n < 00. Finally notice 



where C is some positive constant depending only on Oq. 
Proof: We can decompose Oq as the Fourier series 

0o{xi,X2)= ^ Cj^nej{xi)V2cos{n7TX2) 

j,n>0 

where the convergence holds in L^(0, 1) and 

Y^ \Cj,n\'^ < 00. 
j,n>0 

The solution of ([l])-([5]) starting from Oq then reads 

0{t, Xi,X2)= ^ Cj^rie~^^'^''^''''^*ej{xi)V2cos{n7rx2) 

j,n>0 

and in particular 

0r{xi,X2) = ^ Cj^ne~''^'^'^ "" ^^6^ (xi) cos(n7rx2). 

j,n>0 

This series obviously uniformly converges since 



ari:©'^-^ 



Cj^ne~^^'^''^''^^^ej{xi)V2cos{n7rx2) 



<Cie- 



Moreover 



if we set 



Ro:= < 



4 ^^ 2R -^' 

^ if ^ < 1 
2 ^^ 2R ^ ^' 



-(xi,X2) = V2Y, c,,,e-(^^+-'-')-e,(xi) 

j,n>0 



E(-ir 



i>0 



{2i)\ 



Collecting the three previous estimates yield 

mP (20!i n!f ^ , 



= ^e-^^-e,(xi)^ 



y,2i 



j>0 



i>0 



(2i)! 



2i>ni>0 



i?^ i?f i?2 



j>0 2i>n 



v^(-i)'E^^-'-^ 



-n^TrV (^^) 



2i 



n>0 



(2i)! 



As the two series in the rhs converge the claim is proved. 



■yj,i 



9.(a;) = V2^c„e-"'-'^^(-l) 



n>0 
^2i 



i>0 






i>0 



X 

(2i)! 



= E7i7vh^(-i)'E-»^-"'^^(-)^ 



2>0 



steers the system from the initial state Oq at time to the 
final state at time T. 

Moreover u is Gevrey of order s in t and 1/2 in xi 
on [0,T] X [0,L]; the solution of ([T])-([5]) belongs to 
C{[0,T],L'^{n)) and is Gevrey of order s in t, 1/2 in xi 
and 5/2 in ^2 on [e, T] x H for all e G (0, T). 

Proof: By lemma [T] the state reached at time r is 



The change in the order of summation will be justified once 
we have proved that ^j,i, ^ > 0, is absolutely convergent and 



^\yj,i\ 



i>0 



y.2i 



(2i)! 



< OO, \/X2 > 0. 



For i > let hi{x) := e'^'^^^^Trx)'^' and Ni := (:^) - 
The map hi is increasing on [O, {:;^) ^] and decreasing on 

[(^)' ,+oo) hence 

Y^ hi{n) < / hi{x)dx + hi{Ni) 



n>0 



+ /ii(7Vi + l)+ / hi{x)dx 



/>oo 






hi{x)dx 



<C2 



II 



hi{x)dx; 



r'^Vi Jo 

C2 is some positive constant and we have used Stirling's 
formula. On the other hand integrating by parts yields 

/ hi{x)dx = / hi-i{x)dx 

Jo 2r Jo 






<C3—y=, 

where C3 is some positive constant and we have again used 
Stirling's formula. As a consequence 

|y,,,| < V2 sup \c,^n\Y.K(n) < C(l + ^)4 
where C is some positive constant. Finally 



i>0 



>0 



where the series in the rhs is obviously convergent. ■ 

Theorem 2: Consider any So G L'^(Vt), T > 0, r e (0,T) 
and 5 G (1,2). Then there exists a sequence {yj)j>o of 
Gevrey functions of order s on [r, T] such that the control 



u{t^xi) := < 







A... .. ^ ^?W 



ij>0 



^^-^^^(2^-1)! 



if < t < r, 



if r <t<T 



0r{xi,X2) = ^6 ^'^ej{xi)Yyj, 
i>o 



y,2i 



i>0 



(2i)!' 



with 



%,, = V2^C,•ne-"'-'^(- 



2_2^^ 



-n TT 



n>0 



Clearly the functions 

yj{t):=V2Y,Cj 



i>0 



rrC^) 



satisfy y, {r) = yj^i for all i^j > 0. Hence the desired 
functions yj are given by 



yjit):=<Ps{^,^)y,{t) 



((j)s is defined at the end of the introduction). Indeed by 
lemma [Tithe growth property of the yj/s implies the y^'s 
are analytic on \r^T] hence Gevrey of order s > 1; since 
(j){t) := (j)s ( 1^ ) is also Gevrey of order s on [r, T], so is 
the product yj of ^^ by (j). 

On the other hand the definition of ^^ implies the yj's 
satisfy for alH, j > 



yfir) 
yf(T) 



yj,i 
0, 



SO that 6>(r+, xi,X2) = Or {xi^X2). 

The proposed control u achieves the steering to zero: the 
solution 6 of ([l])-([5]) is then given by ^ and obviously 
satisfies 0{T,x 1^X2) = 0; by proposition [T] it is Gevrey 
of order s in t, 1/2 in xi, and s/2 in X2 on [r^T] x r^. 
Moreover u is also Gevrey of order 5 in t and 1/2 in xi 
on [0,T] X [0,L]. 

It can be shown that for e < t < r the solution 6> of ([T])- 
([5]) is Gevrey of order 1 in t, 1/2 in xi, and 5/2 in X2 
on [r, T] X 1]. The proof is omitted since this fact is not 
used to design the control. To complete the proof we then 
only need to check the compatibility at t = r, ie 

d^0{r+,Xi,X2) = d^0{r,x,,X2) 



for all /c > and x e Q. Indeed using the notations in the 
proof of lemma [T] 

d^0{r^,xuX2) 



E^.(-i)7Sv^*[^"^M'(^ 



i,j>0 



m 



t = T 



i,j>0 ^ ^ m=0 ^ 

E ^^■(^i)(% E (^^ 







(-A,)™e-^^-%,i+fe_™ 



(-A,)™e- 



i,j>0 

x\/2 ^ 9,„e-"''^'^(-n27r2)*+'=-"' 

n>0 
j,n>0 "'^^ ^ ^* 



2i 



i>0 



E 



2_2\fe- 



(-A,)™(-n^7r^) 



n^7r^)^Co ^76" 



-(Aj+n^7r^)r 



= E (-^. 

j,n>0 

xej(xi)A/2cos(n7rx2) 

Hence the solution of ([l])-([5]) belongs to 
C{[0,T],L^{n)) and is Gevrey of order s in t, 1/2 
in xi and s/2 in X2 on [s,T] x 1] for all e G (0,T). As a 
consequence u is Gevrey of order 5 in t and 1/2 in xi on 

[0,T]x[0,L]. ■ 

IV. Numerical experiments 

We illustrate the approach on a numerical example. The 
parameters are L = 1, T = 0.3, r = 0.05, 5 = 1.65. The 
initial condition is the "double step" 



6^0(^1,^2) := < 



f-1 if (xi,X2)G(0,^)x(0,^), 
1 if (xi,X2)G(0,^)x(^,l), 
1 if (xi,X2)G(^,l)x(0,^), 

[-1 if (X1,X2)G(^,1)X X(i,l). 

Its nonzero Fourier coefficients are 



C2^+l,2p+l 



l,p>0. 



^2(2/ + l)(2p + l)' 

Notice 60 is not continuous and that its Fourier coefficients 
decay fairly slowly. The series ([6])-(|7]) and Fourier expansions 
have been truncated at a "large enough" order for a good 
accuracy. 

The following figures are snapshots of the time evolution 
of the temperature from t = Otot = T— 0.025 by increments 
of 0.025 (the snapshot at t = T, being exactly 0, is omitted), 
followed by the complete time evolution of the control effort. 
An interesting question to study in the future is the tradeoff 
between r and T — r to reach the zero state with the smallest 
control effort: longer regularization r or longer duration T—r 
of the active control? 





# 



Control effort u 






^^^ 
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